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ABSTRACT 

We study the structure of shocks in clumpy media, using a muhifluid formahsm. As 
expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous sub- 
shock appears in the structure. This has significant impHcations for the survival of 
dense clouds in regions overrun by shocks in a wide range of astrophysical circum- 
stances, from planetary nebulae to the nuclei of starburst galaxies. 

Key words: Hydrodynamics - Shock waves - ISM: clouds - ISM: globules - ISM: 
kinematics and dynamics 



1 INTRODUCTION 

Astrophysical media are often strongly clumped, with dense 
clouds or bullets of cool gas distributed in a far hotter ex- 
ternal medium. The flows in these media are complex, par- 
ticularly when they are perturbed, for example, by outflows 
from nearby young stars or supernovae. 

Insight into these structures be gained by approximate 
treatments, including multiphase hydrodynamics. Many 
mass-, momentum- and energy-trans fer processes between 
phas e s will be import ant in real flows (Shu et al. 1972, Drew 



of the clouds or as a redistribution of their mass spectrum 
depends on details of radiative and transport processes not 
included in these simulations. 

If the shocks are broadened substantially, the imping- 
ing hot material may be treated as a gradually accelerating 
stream. At least in radiative gas, steady stream interactions 
may be less destructive than those with rapidly accelerating 
shocks (Schiano et al. 1995), as a result of sacrificial strip- 
ping of small instabilities and of the dissipation of frictional 
heat input. Indeed, shocks may be able to cause the aggl om- 



1983; Kamaya 1997b); similar treatmentriire'^seTToirinuF 
t iphase '^pac c plasma flows, as recently reviewed by Szego 
et al. (2000). The simplest cases of pure mass-loading or 



eration of th e condensed material into larger clumps (Ka- 



maya 1997b). The increased pressure of the diffuse gas can 



momentum transfer are interesting limits. 

The mass-loading approximation treats the dense cool 
phase as a distributed source of matter within a hot matrix. 
It is assumed in this treatment that subsequent to the mass- 
loading, dissipative processes such as thermal conduction. 



also lead to condensation of material into the dense phase, 
although this depends on details of the thermal evolution. 

The results of an interaction with many sources has in 
general been treated in a continuum approximation: simu- 
lations in which the interaction of the flow over numerous 
obstacles is studied in detail are only just becoming feasible 



enhanc ;d by turbulent mixing, mean that the hot phase can 
be treated as a single, smooth tiow. I'his approximation has 
been applied, by ourselves and others, to diverse astrophysi- 
cal objects, such as Wolf-Rayet and planetary nebulae, ultra- 
compa ct Hii regions, starburs t galaxies and active galactic 
nuclei (Chevalier fc Clegg 1985 ; Hartquist et al. 1986t Arthur 



.1^ 



un, Jones & Norman 1996; Poludnenko, Frank & Black 



et al. 



1997 



1994 



Redman et al. 1996 ; Smith 1996 ; Wilhams et al 



Underlying these approximations are the details of the 
interaction between a dense clump and a steady wind or a 
shock have been studied numerically by many authors (e.g. 



Woodward 1972; Nittman, Falle & GaskeU 1982 



Klein, Mc- 



Kee fc Colella 19941 ). Where a strong shock propagates over 



a dense cloud, these authors find that the cloud is crushed 
in the direction parallel to the incoming shock, and ex- 
pands in the lateral direction. Whether the net effect on 
the cold phase gas can better be described as destruction 



2QQ1D , and the number of obstacles in these studies 
is limited. In Fig. 1 of Jun et al., the leading shock is ap- 
parently slowed and weakened somewhat by the presence of 
the clouds, but it is difhcult to infer the overall behaviour 
of the shocks from this paper as only two simulations are 
illustrated. The finite resolution of such studies leads to nu- 
merical viscosity and numerical thermal conduction which 
can result in far more effective coupling between the phases 
than might be expected in reality. 

In the present paper, we treat the multiphase structure 
using a momentum-loading formalism. This allows a con- 
tinuum treatment of the internal structure to maintain the 
distinction between cool- and hot-phase properties, while in- 
cluding both as dynamically active components of the flow. 
We parameterize the uncertainty about the detail of the 
momentum-transfer process by considering separate forms 
of the coupling terms, and look for general features of the 
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results. We study the structure of steady shocks in a multi- 
phase ISM in this approximation, for the case in which the 
two phases ahead of the shock are comoving. This case is 
somewhat distinct from the formation of termination shocks 



in an inpinging wind driven by mass injection from gas 



which has been overrun, e.g. in syst ems such as cometary 
particle injection in th e solar wind (Biermann et al. 1967 
Zank & Oughton 1991), or clumps which have been overrun 



by bubbles dri ven by stellar wind s ( Williams et al. 1995 ) 



or supernovae (Cowie et al. 1981 ; Wolff fc Durisen 1987 
Pittard et al. 200l| ). 



Ka- 



K amaya has studied the wave dispersion relation 
maya 1 )97ar) and time-dependent shock tube problem (Ka- 
maya 1 )97b[ ) for two phases initially at rest with respect to 
each other. The time dependent solutions in the latter paper 
were not integrated for a time sufficient for the two phases 
to reach equilibrium between the shock and the rarefaction; 
only one case (with a leading subshock) is presented. Our 
study complements this work, by determining steady-state 
structures for a well-defined, if simplified, set of model equa- 
tions. 

We discuss the implications of our results for various 
astrophysical systems. 



2 BASIC EQUATIONS 

We consider a system of two phases, a distributed adiabatic 
hot phase and a clumped cold phase which fills a negligi- 
ble fraction of the volume of the fiow. In reality there will 
be significant substructure in both the hot and cold phases, 
but we use continuum equations for both phases on the as- 
sumption that a suitable coarse-graining process can provide 
consistently-defined local mean values for the flow density, 
velocity and pressure. 

We assume that the hot phase obeys mass, momentum 
and energy equations 



_d 

da; 



_d_ 

da; 



IP 

7-] 







= <S [t;i + a{v2 — vi) 



(1) 
(2) 

(3) 



where pi is the gas density in the hot phase, m is the flow 
velocity, and p is the pressure. The source term S is the slip 
force between the two phases, and a is the fraction of the 
frictional energy dissipation resulting from this slip which 
acts to heat the hot phase. In the absence of a second phase 
(iS = 0) , these equations are simply the equations of steady- 
state flow for a fluid with adiabatic index 7. 

The clumped cold phase obeys the mass and (ballistic) 
momentum equations 







d^ ^P'""') 



= -<s. 



(4) 
(5) 



where p2 is the mean density per unit volume of the cold 
gas, rather than the density within individual clumps. These 
equations form an autonomous system, as x appears only in 
the differential operators. 



A particularly simple set of equations is used here, con- 
sidering only momentum-transfer effect. To accurately treat 
a specific flow, a wide range of ma ss and energy-transfer 



processes should also be considered (Shu et al. 1972; Wolfi 



& Durisen 1987 ) . However, the treatment of the astrophysi- 



cal plasma as consisting of locally uniform interacting media 
is itself a considerable approximation, and we will see that 
our further simplification of the system leads to interesting 
results. 

Equations (^) , (^) and the sum of equations (^ and (^ 
can be immediately integrated, to give 



$1 = pivi 
$2 = P2V2 

P + P1V1+P2V2 = n. 



(6) 
(7) 
(8) 



In addition, <S can be eliminated from equations (|l|)-(^ to 
give 



d 

da; 



-ypvi 



L " / 3 , 3\ 
J + 2 "1 + P2V2) 



(1 — a)vi 



dp 
dx 



(9) 



From this equation it is clear that for a = 1, the total energy 
fiux 

£^^ + l{piv',+P2vl) (10) 

7 — i Z 

is constant. For a = 0, the hot fiow is adiabatic so long as 
there are no discontinuities in the flow, and hence p oc v^^ . 

For general a, the dependence on x can easily be elimi- 
nated, as can those on pi, p2 and p, using the integrals (|^)- 
(^). If the solution is smooth, the result is a first-order o.d.e. 



dwi _ Vi — 0(7 — 1)(«2 — Vl) 

dv2 7(II/<I>2) - (7 + l)(<E>i/<&2)wi - 7^2 ' 



(11) 



This equation is integrable for all a. However, the cases a = 
and Q = 1 are sufficient for our purposes here. We note 
that this equation depends only on the fraction of the energy 
deposited in each of the phases, and not on the form of the 
momentum coupling. 

By including the momentum coupling, the full shock 
structure as a function of x can be derived for each solution 
of equation (pT|), using 

Since <S is a function only of local conditions in the flow, it is 
clear that V2 must decrease monotonically through the shock 
structure (unless S — somewhere within the structure). 

Dividing equation ( [ll| ) by equation (|l^, we flnd that 
a critical point exists, as expected, when the flow is at the 
hot-phase sound speed, i.e., vf = jp/pi, or 



Vl = 



7 n — ^2V2 

j+1 ^1 



(13) 



It may be possible for a continuous solution to pass through 
this critical point if <S = 0, or if 



■y2 = 1 + 



0(7 - 1) 



Vl 



(14) 



(for a / 0). 
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2.1 Shock solutions 

Shock solutions for the basic equations exist where the flow 
velocity into the discontinuity exceeds the sound speed in 
the well-coupled limit, i.e. 



Vi > c 



7Pi 



Pl,i + P2,i 



1/2 



(15) 



where pi,i and p2,i are the upstream densities, pi is the up- 
stream pressure, and Vi — vi^i = U2,i is the upstream velocity. 

For c* < Vi < ci.i = {"fpi/ pi,i)^^^ , the shock structures 
will be continuous in both fluids (C-type). For higher veloc- 
ities the shocks will be led by viscous subshocks in the hot 
gas (J-type) (following the nomen clature defiri ed for shock 
waves with magnetic precursors by |Z)raine 198C ) : as the hot- 
phase sound speed is the fastest characteristic speed in the 
undisturbed gas, the subshock must occur at the start of the 
resolved shock structure. For J-type shocks, the velocity of 
the hot phase just after the subshock is 



7 - 



1 27 Pi 

-Vi + 



(16) 

the standard 



7 + 1"' ■ 7+l<&i' 

and the pressure is ps ~ Pi + 4>i(fi 
Rankine-Hugoniot conditions. 

For a — 1, energy conservation gives that the final ve- 
locity, V{, of a shock with upstream velocity Vi and pressure 
Pi is given by 

27 



Vi 



7 ■ 



-Vi + 



Pi 



(17) 



7 + 1 ■ 7 + 1 $1 -f- $2 

whether or not it contains a viscous subshock. 

For Q = 0, the situation is more complex. For C-type 
shocks, V[ is given by 



Oi (^^j + ($1 + $2)l'f = Pi + ($1 + $2)"i. 



(18) 



For J-type shocks with Vi > ci,i, however, the final velocity, 
Vf, is given in terms of v^, Ps behind the subshock by 

Ps f — ) + ($1 + $2)l'f = Ps + $l«s + 4'2Wi. (19) 



The expressions for V{ are different for C- and J-type shocks 
for all a 7^ 1. 

Note that no shock solution can reach a sonic point 
where equations (|l^) and (|lj) are satisfied. We can see this 
because in order to satisfy these conditions. 



V2 = V2 



[pi + ($1 + <i>2)vi] I 

1 + 7\ / a(7 



l + a(7- 1) 



+ $2 



(20) 



If 7 > 1 and < a < 1, > ui, but 112 must decrease 
monotonically through the shock structure [see the discus- 
sion after equation (nd) above] so no such critical point can 
occur in practice. 

2.2 Drag force 

The spatial structure of the shocks which we calculate is de- 
termined by the drag force, S. The simplest assur nption is a 
relat ion which takes the form of Stokes' law drag (Batchelor 
1967|), for which 



XpxP'zivi - V\). 



(21) 



If the cool phase consisted of isolated, rigid, spherical clouds, 
then A = (S-kR^v jM^, where 7?c and Mc are the radius and 
mass of the clouds and v is the kinematic viscosity of the in- 
tercloud medium. Interpreted in this direct sense, this law is 
only valid where the Reynolds numbers 27?c|«2 — v\\lv "^-y 1, 
which will not apply in astrophysical situations: at larger 
Reynolds numbers, the effective drag coefficient increases. 
While it remains of interest to consider this Stokes' drag 
force for an arbitrary drag coefficient, the second order term 
in the expansion for the drag force will cancel the depen- 
dence on fiow viscosity but introduce a higher order depen- 
dence on the slip velocity. 

Assuming that the second-order term dominates leads 
to a drag law proportional to th e square of the fiow veloc- 
ity, as used by Shu et al. (1972), a dependence which also 
corresponds to the 'ballistic' case, where the drag force is 
calculated by assuming independent collisions of hot-phase 
particles with the dense clouds. It agrees well with empirical 
data for rigid projectiles. In the limit which we are consid- 
ering, of small filling factor of clouds, the drag force may be 
written 



S = C 



97r72 
128 



1/3 



p2(v2 — Wl) 



Pi /^£2\ 



4/3 



\V2-Vi\, (22) 



where the typical mass of a cold clump is Mc and the typical 
internal isothermal sound speed is C2. As above, ci is the 
adiabatic sound speed in the diffuse medium. 

An alternative form for the drag force is based on the 
work of Hartquist et al. (1986, hereafter HDPS), who took 
into account the effects of deformation of the clouds. These 
authors considered the dynamics of a fiow around a clump in 
both the subsonic and supersonic limits, assuming that the 
interaction zones were small enough that each clump inter- 
acted with a smooth upstream fiow. They determined forms 
for the mass-loss rate of the clumps, which have been used 
in subsequent papers to predict a smoothed mass-input rate 
for the global fi ow. a nd agree well with the numerical results 
of Klein et al. ( 1994 ) for strong shocks propagating over iso- 
lated clouds. It is a simple matter, however, to reinterpret 
their results as a momentum transfer: the value is just the 
mass loss rate between the phases derived by HDPS mul- 
tiplied by the slip velocity. If the slip speed is slower than 
the sound speed in the hot phase, we find that the coupling 
force is given by 



P2{V2 - V\) 



P\\V2 — V\\ci 



1/3 



\{\v2~vr\lc^f^\\) 



We see that in both the Shu et al. and HDPS formulae, the 
momentum transfer rate is proportional to the mass density 
in the cool phase divided by the clump mass to the power 
1/3: if the cold gas is placed in larger clumps, the ensemble 
is more permeable to the hot phase fiow. 

It is of interest to study the results of these forms of 
slip term in the test-particle limit (i.e. to treat the case in 
which phase 2 corresponds to a single cloud, and the hot 
phase velocity is effectively constant). From equations (^ 
and (jsj), we find that the velocity of the cloud relative to the 
hot phase gas is 



di;2 
"dF 



$2 



(24) 
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(a) 



(b) 



a = 0; = 1 



a = 0; = 10 





(c) 



(d) 



a = 1; = 1 



a = 1; = 10 





Figure 1. Variation of vi and V2 tlirough resolved shocks for a = 0, 1 and density ratios <I>2/'I'l = 1, 10. Structures for shock Mach 
numbers vi/ci j = 2, 1.2,0.8 are shown on each plot. Circles show the initial and final states and the solid curve the resolved structure, 
while the arrow shows the viscous subshock in the structure if one is present. The dotted curve shows the equilibrium condition vi = V2, 
while the dashed lines show the critical conditions vi = ci (for each value of fi/ci;) and, for a = 1, equation (^). 



For Stokes' drag, fH = 1 and A = Api, while for the (sub- 
sonic) HDPS law, (3 — 8/3 and A is a constant independent 
of V2. Integrating, we find that for (3 — 1, 



V2 = Viexp{—At) 



V2 = 



Vi 



-^[1- exp(-At)] , 



for /3 = 2, 



Vi 



1 + Avit 



-log(l + A^;it), 



(25) 
(26) 



(27) 
(28) 



In the Stokes' drag case, the dense cloud tends to a finite 
displacement, Vi/Xpi, from its initial position in the diffuse 
flow. However, for the subsonic HDPS law (and in general 
for all /3 > 2) there is no such limit. 

This observation relates to d'Alembert's paradox - the 
prediction of zero drag on a body in irrotational flow, con- 
trasted to practical experience in which turbulence in the 
wake of bluff bodies maintains drag. For slender bodies in 
flows of hi gh Reynolds nu mber, very low drag is observed 
in reality (Batcholor 1967), and so long slip distances could 
indeed occur for in some cases. We include results for both 
types of behaviour in what follows. 



and otherwise 



V2 = Vi[l + {f3~l)Av';'''t 



-l/(/3-l) 



(/3 - 2)A 



{[1 +{p-i)Avrhf 



(29) 

- 1 jso) 



3 SHOCK STRUCTURES 

In Figure |l[ we plot the structure of the shocks in drag- 
coupled two phase media in the («i, W2) plane, for a = and 
a = 1 respectively. These plots are independent of the form 
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(a) 



(b) 



= 0; v-/c,- = 0,6; *g/*| = 1 



a = 0; v./c.. = 0,8; ^^''^i = ^'•^ 




10 




(c) 



(d) 



cx = 0; v./c, i = 1.2; = 10 




5 8 




(e) 



(f) 



: = 0; vjc, , = 2; »,/», = 1 



i = 0; v./c, , = 3; = 10 



10 



B 10 



Figure 2. Shock structures for a = and various shock Mach numbers vi/ci^i and density ratios <l?2/"I'i and with a sHp force given by 
equation (|2l|), for ci^i = pi^j = 1, The dotted curves show the initial and final velocities of the shocks, as well as tii for a limiting case 
with V2 = which relaxes to the same final state. The solid curve between the upper dashed lines shows the variation of vi through the 
shock (including the initial subshock for di ; > ci ;), while the dashed curve shows V2- Additional solid curves show solutions relaxing to 
the pre-shock and post-shock equilibrium states. 
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10 



(c) 



(d) 



<* = 1: i = 1.2; = 1 



= 1: Vj/Cii = 1.2; fg,/*! = 10 



Figure 3. Shock structures for Vi/c-i^^ = 1.2 as in Figure Q, except that a = 1, and for (c), (d) the HDPS sUp law is used. 



of the drag force, S [cf. equation (|ll|)]. The two columns 
of these Figures are for a density ratio of pi,i/p2,i of 1 : 1 
and 1 : 10, respectively, with curves for several shock Mach 
numbers (with respect to the adiabatic sound speed in the 
smooth phase, ci,i) shown on each plot. For comparison, 
for a 1 : 1 density ratio, the equilibrium critical speed is 
while for a 1 : 10 ratio is is ci,i/\l^''^. For a = 0, 
the exhaust velocity of the shocks decreases for the strongest 
shocks, while for a = 1, the exhaust velocity is increased as 
a result of the energy deposition in the hot phase. 

For the weakest shocks, we see that the velocities in the 
two phases remain very similar, and have continuous (C- 
type) profiles. For case "l>2/$i = 10, substantial compres- 
sions of the medium can be achieved in these C-type shocks. 
As the shocks strengthen, the smooth phase begins to de- 
crease in velocity significantly before the clumped phase. 

Once > ci,i, a viscous subshock is required in the 
solutions (shown as an arrow), which are hence J-type. It 
is apparent that, for q = 0, increasing clumped-phase den- 
sity results in a stronger shock (for fixed fi/ci^i), and more 
post-shock velocity structure in the smooth phase gas. These 
effects are less marked for a = 1. For shocks with a = 1 



stronger than those shown here, the smooth-phase gas ac- 
celerates after passing through the subshock. 

In Figures ^ and ^ we show the internal structure of the 
shocks as a function of distance, for the Stokes' drag law, 
equation (|2l[), and for the HDPS slip law, equation (p3|). The 
shock structures get narrower as the shock velocity increases, 
while going from q = to a = 1 broadens them somewhat. 
Increased mass flux in the cold phase, $2, appears to narrow 
the shocks. This is a result of the reduction of the critical 
speed for the well-coupled limit, which means that for fixed 
velocity through the upstream gas, the shocks with higher 
$2 are stronger. Note in particular the extremely small slip 
velocity between the phases throughout the C-type shock 
shown in Figure ^a). 

The dashed lines in the Figures show the variation of 
V2 through the shock: it is apparent that the velocity dif- 
ference between the phases increases for stronger shocks. 
Additional solid lines show supersonic and subsonic equi- 
libria. For shock speeds smaller than the hot-phase sound 
speed, the solutions of the o.d.e. diverge from the higher- 
velocity eguihbriurn^ we have previously dis- 
cussed ( Williams fc Dyson 199(j ), the higher- velocity equi- 
librium solution is not physically unstable, but this solu- 
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tion topology will occur if fully resolved shock structures 
are present. 

Comparing the results for Stokes' and HDPS drag, we 
have set the constants of proportionality so the initial post- 
shock relaxation occurs at a similar rate. It is apparent that 
the HDPS form leads to an extended post-shock relaxation 
region, as the small slip velocity in this region leads to de- 
creased effective friction. 



3.1 Vorticity 

The question of whether the passage of a shock through a 
clumpy medium smooths out the clumpy structure depends 
largely on transport processes such as viscosity and ther- 
mal conduction (although in astrophysical plasmas, radia- 
tive processes can also play an important role) . The vorticity 
produced by the passage of the shock is crucial in enhancing 
transport processes, as a result of the formation of viscous 
boundary layers. 

If there is a viscous subshock, then the flow is not ho- 
moentropic and so Kelvin's circulation theorem (e.g. Batch- 
elor 1967, p273-277) does not hold. Shocks in the flow will 
generate vorticity when their strength changes across their 
surfaces. Considering the small scale structures of the flow 
as the subshock propagates through it, it is apparent that 
vorticity will be generated in the hot phase where the inci- 
dent shock refracts around obstacles, and by reflected shocks 
ahead of them. The forward shock driven into each clump 
will also lead to creation of vorticity within the clump. This 
wide distribution of vorticity within the flow will have a sig- 



Ax = 



^^f (i 



(32) 



Jun, Jones & Norman 1996) 



nificant disruptive effect (cf . Klein, McKee & Colella 1994 



However, if there is no viscous subshock and a — 0, the 
hot-phase flow is homoentropic and (as a result of Kelvin's 
circulation theorem), the vorticity produced will be entirely 
confined to the cold phase and the contact between the 
phases. Even for a ^ 0, the generation of vorticity should 
be confined to the regions of the hot phase flow close to the 
boundary layers in which the frictional heating is deposited. 

As the slip velocity gets smaller, the shear velocities 
over the surface of individual clumps decrease, and so does 
the effectiveness with which cool material is entrained from 
their surfaces. Instead, the cool clumps are rather gently 
accelerated as they cross the broad shock structure. 

We can develop various approximate measures of the 
destructiveness of the boundary layers induced by passage 
through a shock. For instance, the amount by which two 
initially coincident particles, one in each phase, are dragged 
apart by the front, 



Vf 



1 

Vl 



-dx. 



V2 



(31) 



where Vf is the final equilibrium velocity behind the shock, 
and Vl and V2 are the velocities of the different phases as they 
pass through it. The value of Aa; may be compared to other 
characteristic lengthscales in the problem. If it is smaller 
than the size of an individual cold clump, then the passage 
of the shock will have very little effect on the structure of the 
medium, while if it is larger than the inter-clump separation, 
the cool phase will likely be dispersed into a fine aerosol. 

Using equation (|) and the Stokes' drag force, equa- 
tion (|l|), we find that 



It is clear that as the shocks weaken or the coupling im- 
proves, Aa; becomes smaller. By contrast, the slip law of 
HDPS leads to increasingly poor coupling between the 
phases at small slip velocities compared to Stokes' drag and 
the slip distance would be expected t o increase for small 
velocity jumps (cf. Yuan & Wang 1982). 

An alternative measure of destructiveness is the max- 
imum slip velocity between the phases. This can be com- 
pared, for instance, to the peculiar velocities of the clumps 
in the upstream multiphase gas: if the maximum slip veloc- 
ity is smaller than this value, it is difficult to see why the 
multiphase structure should be substantially effected by the 
passage of the shock. The interstellar medium has significant 
structure on large scales, while for a purely hydrodynamic 
flow in radiative equilibrium, it would be expected that dif- 
ferent phases would be dispersed as an aerosol. However, 
scale lengths are introduced by effects of such as self-g ravity 
(the Jeans length), thermal conduction [the Field (1961:) 
length] or photoionization (the recombination length), re- 
sulting in large-scale clumping. If these structures are to 
form, then they must have a binding energy sufficient to 
outweigh the higher entropy of the dispersed state. 

It is clear, from Figure |l|, that the maximum slip veloc- 
ity decreases rapidly with shock strength. For example, for 
Q = 0, the maximum value of the slip velocity, V = V2 — vi 
is given by 



V* 



$1 -I- $2 
7$2 



+ 7— -(7+1) 

Vl 



(33) 



for a shock with exhaust velocity velocity V{, and velocity vl 
at the maximum of V* . As — > Wf, V* = 0[{vf — wj')^], so 
for weak shocks the maximum slip velocity decreases rapidly. 

If the square of the maximum slip velocity is less than 
the effective binding energy per unit mass of these struc- 
tures, then the shock should be able to change the velocity 
of the flow as a whole without greatly effecting its multi- 
phase structure. While some details of the shock structures 
vary substantially between different forms of the momentum 
coupling law, it is in all cases that the maximum slip veloc- 
ity decreases rapidly as the overall shock velocities decreases 
below the hot phase critical speed. 



4 APPLICATIONS 

The relatively benign environment of resolved shocks may 
have important implications for the evolution of multiphase 
structures in a wide range of astrophysical circumstances. 
In the present section, we outline some of those we expect, 
pending future more detailed investigations. 

In Figure ^ we show a schematic picture of the struc- 
ture of shock-tube problem, in which two gas components 
collide. The component moving in from the right is smooth, 
while the component moving in from the left has a clumped 
structure. If the flows from either side are sufficiently rapid, 
the solution will take up the usual form of reverse shock, con- 
tact discontinuity, forward shock. In Figure ^ we illustrate 
the situation in which the reverse shock into the multiphase 
medium is of J type: after they pass through the shock, the 



© 0000 RAS, MNRAS 000, 



8 R.J.R. Williams & J.E. Dyson 



Contact 



Shock Discontinuity 








°o 

OqO 




.On-- 




00° 

o^op 




000 

« 




QO 
^ 






0^ 


00 
oo ° 





Sliock 



Figure 4. Schematic drawing of the structure of a multiphase 
shock. 



clumps become elongated and are eventually destroyed (al- 
though the deformation of the clumps will obviously be far 
more complex than that shown in the Figure). If the shock is 
of C type, however, we have argued that the clumps will be 
far more gradually decelerated and as a result could remain 
as fast-moving coherent structures within the protective co- 
coon of the shocked gas layers. 

If a clump slips forwards through the contact discon- 
tinuity, its environment will become somewhat more harsh. 
This will occur if the dense clumps slip through the hot 
material for a distance larger than that between the shock 
and the contact discontinuity, or, in an environment which 
is not plane-parallel, as a result of sideways escape of the 
diffuse gas. As the forward shock will often be curved, the 
flow beyond the contact discontinuity will often have signif- 
icant shear flows. Once a clump passes through the line of 
the forward shock, however, it will shock directly against the 
supersonic impinging material and will rapidly dissipate. 

This may help expla in the survival of knots and bu llets 
in regions su ch as Orion ( Tedds^Jgrand fc Burton 1999 ) and 
Eta Carinae (Meaburn et al. 1996). In systems such as these, 
the spatial density of the bullets may be high enough that 
they maintain a relatively shared bow-shock, and so reduce 
the destructive vorticity of the post-shock flow. 

It is generally thought that the structures of planetary 
nebulae are formed when a hot, diffuse wind from a white 
dwarf remnant overtakes a cool, dense stellar wind blown 
off during the stars previous evol u tion along the asymptotic 



giant branch (Kwok et al. 1978; Mellema 1995). The cool 



such as the Abell 
start to interact with the hot stellar wind 



wind is likely to be highly structured, and clumps advected 
into t he hot bubble can have important effects on the global 
flow ( Arthur et al. 1994|), as well as being obser ved in objects 
30 ( porkowski et al. 1995| ). Before they 

these clumps 

will have had to pass through an outer shock without being 
destroyed. In particular, the outer region of shocked gas may 
be broadened substantially, compared to numerical models 
which have often assumed the cool wind is uniform. 

This analysis also provides insight into the structure of 
the flows around nuclear starbursts (Strickland & Stevens 
2000, and references therein). Energy input from supernovae 
and stellar winds in these starbursts drives plumes of hot 
gas perpendicular to the plane of the host galaxy, over dis- 
tances of many kiloparsecs. In the plane of the galaxy, the 



hot wind from the starburst interacts with the ISM. To date, 
simulations have in general assumed that the host galaxy 
ISM may be treated as a smooth gas, at a temperature of 
~ 10^ K. This temperature is similar to the mass weighted 
mean random velocity in the clumpy ISM, and so gives a 
fairly reasonable impression of the coarse-grained structure 
of the ISM before the intera ction with the nuclear wind be- 
gins. Strickland & Stevens (2000) do investigate the effects 



of mass-loading in some of their models, and flnd relatively 
small effects on the flow structure. However, in the galactic 
plane in their simulations the distributed gas has a temper- 
ature of only ~ 7x 10* K, while we flnd here that the hotter 
components of the gas will be very important for transmit- 
ting pressure waves. 

When the nuclear starburst starts to drive a flow out- 
wards in the simulations, the region around the starburst is 
first evacuated of the ambient ISM gas. After this, a nuclear 
wind develops which becomes supersonic close to the edge 
of the starburst mass injection. The resulting nuclear wind 
shocks against the residual ISM with the usual structure of a 
forward-shock, a contact discontinuity and a reverse shock. 
Some cool gas may be dragged away from the galactic plane 
by shear flows along the contact discontinuity, but in general 
the structures seen in the galaxy plane are very sharp and 
well defined. 

In contrast, the broadening of the shocks will mean that 
clouds in the ISM will be rather gradually driven away from 
the starburst by a diffuse, broad outer shock. The small force 
on the clouds will result in a secular increase in orbit energy, 
rather than any dramatic outflow. As a result, in the plane 
of the galaxy the inner shock and sonic transition around the 
starburst may take a considerable period to develop. The hot 
wind from the central nuclear starburst will be surrounded 
by a weaker wind which has percolated through the ISM, and 
which carries some shreds of advected cool material along 
with it. 



5 CONCLUSIONS 

In this paper, we have discussed the effects which the mul- 
tiphase nature of a medium can have on the structure of 
shocks passing through the medium. The shocks rapidly be- 
come broader, particularly when the shock speed in the well- 
coupled limit is smaller than the sound speed in the medium 
which fills most of the space. We discuss how, in this case, 
multiphase structure would be expected to survive the pas- 
sage of the shock. 

Broadened multiphase shocks will be important in a 
wide variety of astrophysical circumstances, as we have 
discussed above. Equations such as our equations (|l|)-(|^, 
which model the effects of small-scale structures, can be used 
in time-dependent hydrodynamical modelling of large-scale 
structures. For a small computational expense, significant 
additional insights into the structure of astrophysical flows 
are available, without the need to resolve the structure of 
the flow to its flnest details. 
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